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Abstract. We prove strong convergence to singular limits for a linearized fully 
inhomogencous Stefan problem subject to surface tension and kinetic under- 
cooling effects. Different combinations of cr — > ctq and S So, where <r, ctq > 
and <5, <5rj > denote surface tension and kinetic undercooling coefficients re- 
spectively, altogether lead to five different types of singular limits. Their strong 
convergence is based on uniform maximal regularity estimates. 

1. Introduction 
The aim of this note is to consider the fully inhomogeneous system 

(d t -cA)v = f mJxi n+1 , 

7« ± — aA x p + 6dtp = g on J x M™, 

dtp + lcd y (v - ap E )] = h on J x R™, (1) 

v(0) = v a in K"+\ 

p(0) = po inM", 

which represents a linear model problem for the two-phase Stefan problem subject 
to surface tension and kinetic undercooling effects. Here 

denotes the temperature in the two bulk phases = {(x,y); x £ R™, ±y > 0}, 

and we have set M" +1 = M" +1 U K" +1 and J = (0, T). The function p appearing in 



1991 Mathematics Subject Classification. Primary: 35R35, 35B65, 80A22; Secondary: 35K20. 
Key words and phrases. Stefan problem, free boundary problem, phase transition, singular 
limits, maximal regularity. 

The second author is supported by the Center of Smart Interfaces at TU Darmstadt. 

1 



2 



JAN PRuSS, JuRGEN SAAL, AND GIERI SIMONETT 



the boundary conditions describes the free interface, which is assumed to be given 
as the graph of p. We also admit the possibility of two different (but constant) 
diffusion coefficients c± in the two bulk phases. The parameters a and 5 are related 
to surface tension and kinetic undercooling. The function pe is an extension of p 
chosen suitably for our purposes. Here it is always determined through 

(d t -cA)p E = in J x R n+1 , 

1P E = P on J x R™, (2) 

p E (0) = e-^-^po inR" +1 . 

Using this notation, let \cd y (v — pe)\ denote the jump of the normal derivatives 
across R™, that is, 

\cd y {v - pe)] ■= c +1 d y {v + - p%) - c-'ydyiv' - p E ), 

where 7 denotes the trace operator. The coefficient a is supposed to be a function 
of 5 and a, that is, a± : [0, oo) 2 — > R, [(8,<r) H> a±(5,<r)]. It is further assumed to 
satisfy the conditions 

o± G C([0,oo) 2 ,R), a±(0,0)>0. (3) 

Recall from [10] that the introduction of the additional term 'ape' with a± > 
in the situation of the classical Stefan problem is motivated by the following two 
facts: for suitably chosen a (depending on the trace of the initial value and d y pE) 
it can be guaranteed that a certain nonlinear term remains small for small times. 
On the other hand, the additional term 'ap^' is exactly the device that renders 
sufficient regularity for the linearized problem. Note that, concerning regularity, 
this additional term is not required if surface tension or kinetic undercooling is 
present. However, in order to obtain convergence in best possible regularity classes 
for the limit tr, 5 — > 0, we keep the term 'aps' in all appearing systems. Since 
the data may (in general even must; see Remark 1) depend on a and 8 as well, 
a is a function of these two parameters. The natural and necessary convergence 
assumption (10) then implies that we can assume that a± € C([0, oo) 2 , R). This 
continuity will be important in deriving maximal regularity estimates for related 
boundary operators; see the proof of Proposition 1. 

The results of this paper on system (1) represent an essential step in the treatment 
of singular limits for the nonlinear Stefan problem on general geometries. This will 
be the topic of a forthcoming paper. 

To formulate our main results, let W* (R ri ), s > 0, p £ (1, 00), denote the Sobolev- 
Slobodeckij spaces, cf. [15] (see also Section 2). Depending on the presence of 
surface tension and/or kinetic undercooling we obtain different regularity classes 
for p, the function describing the evolution of the free interface. To formulate this 
in a precise way we define for J = (0, T) and 5, a > 0, 

E 2 T (S, a) := [p G E|(0, 0) : d|M| E| (1 , 0) + £r||p||^. ( o,i) < 00} , (4) 

equipped with the norm 

II ' HeI^.ct) : = II • IIm|(o,o) + <5|l • IIe|,(i,o) + v\\ ■ IIk|(o,i)> ( 5 ) 
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and where 

i£(0, 0) := W^ 2 - x ' 2p {J, L P {W 1 )) n W X (J, W^-^iW 1 )) n L P (J, V^-^^M")), 
E|(1,0) := W 2 - 1 / 2p {J 1 L p (R n ))nW x (J 1 W 2 - 1 / p (R n )) 1 

E|(0, 1) := W^ 3/2 - 1/2p ( J, L P (K")) n W^ 1/2p ( J, W 2 (R n )) n L p ( J, ^^(R")), 

equipped with their canonical norms. For the different values of S and cr (i.e., 
S = a = 0, or <S > and cr = 0, or 5 = and cr > 0, or 5 and cr > 0) we obtain four 
different regularity classes for p. This leads to the following five types of singular 
limits for problem (1): 

(1) (5,(7)^(0,0), 5,CT>0, 

(2) (S, a) -> (5 , 0), for <5 > fixed, 

(3) (<5, a) -> (0,cr o ), for cr > fixed, 

(4) (*,0)-> (0,0), 

(5) (0,<7)-> (0,0). 

Our main result, Theorem 1.2, covers convergence results for all these limits. 
In the sequel 

f 1, *>0, 
Bg(*):=< 0, i = 0, 
I -1, t<0, 

will denote the sign function. Our first main result is on maximal regularity Here 
we refer to Section 2 for the definition of the space of data Ft {5, cr). The essential 
difference to corresponding results in previous publications is the uniformness of the 
estimates with respect to the parameters S and a. 

Theorem 1.1. Let 3 < p < oo, R,T > 0, < S, a < R, and suppose that 
a = a(S,a) is a function satisfying the conditions in (3). There exists a unique 
solution 

(v,p, PE ) = (v^\p^\p^ } )eE T (S,a) 
for (l)-(2) if and only if the data satisfy 

(f,g,h,v Q ,p ) € ¥ t (S,<t), (6) 

rfi ~ aA xPo + 6 (h(0) - ind y {v Q - ae-l^l^-^^po)]) = g(0), (7) 
and, if 5 = 0, also that 

a(h(0) - [cydyvo]) G W 2 ~V p (R n )- (8) 
Furthermore, the solution satisfies the estimate 

\\(v,p,pE)\\m T (s,a) < C (||(/,5,ft,uo,^o)||F r (o,o) + ( <5 + <T )IIPo|l w 4-3/ P(R „ ) 

+ <r\\h(°) - lcrfd y vol\\ w 2-6/ P{Kn) ^ , (9) 
where the constant C > is independent of (5, a) E [0,R] 2 . 
Our main result on convergence of singular limits is 
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Theorem 1.2. Let 3 < p < oo, R,T > 0, < 5 < 5 < R, < <r < a < R, 
and a = a(8,o~) be a junction satisfying the conditions in (3). Set fi := (5, a), 
fiQ := (<5o,<to), and Iq := [i5o,-R] x [fo,-R]- Suppose that 

and that the compatibility conditions (7) and (8) in Theorem 1.1 are satisfied for 
each [i G Iq. Furthermore, denote by (y 11 , p^,) the solution of (l)-(2) given 
in Theorem 1.1 that corresponds to the parameter /i = (<5, a) G Iq. Under the 
convergence assumptions that 

(r,9^h^v^p^^(f^,g^,h^,v^,p^) m F T ( Mo ), (10) 

and, if So = 0, that 

a (^(o) - [cydyvS]) -+ <jo (h(or - [c 7 a^CD m w^(^ n ) (n) 

and, if S = o-q = and 6 > 0, also that 

(S + a)p^0 in Wp~ 3 / p (K") (12) 
on the data, we obtain strong convergence of the solution, i.e., we have that 

K/,^)^^,^^^) m Et(mo). (13) 

Remark 1. (a) Note that for 8 > condition (8) follows automatically from 
condition (7). 

(b) Conditions (10) and (12) for the last component are obviously satisfied for a 
fixed initial interface in Wp~ 3 ^ p (R™), i.e., if we assume pg = po G Wp~ 3 ^ p (W l ) for 
all fi G Iq. But observe that, due to condition (7), it is not possible to fix vo as well. 

(c) In analogy to (a) note that for 5q > assumption (11) follows automatically 
from (7) and (10). Also observe that in the case S = 5q = ctq = condition (12) 
follows automatically from conditions (7) and (10). 

(d) In the case Sq = a = conditions (11) and (12) express that ||pg|| and 
\\h(0)i*-l'yd y vg}\\ might blow up in Wp~ 3/p {R n ) and Wp~ 6/p (R n ) respectively, but 
slower than a and 6 tend to zero. This seems to be natural in view of the fact that 
we do not have p °' 0) = p (0 ' 0) | t=0 G Wp~ 3/p (W l ) and 

M0)(°<°> - Indyv^j = V°' 0) k=o e w*- 6 / p (R n ) 

from the regularity of solutions in the situation of the classical Stefan problem. 

The Stefan problem is a model for phase transitions in liquid-solid systems that 
has attracted considerable attention over the last decades. We refer to the recent 
publications [5, 10, 11, 13, 12] by the authors, and the references contained therein, 
for more background information on the Stefan problem. 

Previous results concerning singular limits for the Stefan problem with surface 
tension and kinetic undercooling are contained in [1, 16]. Our work extends these 
results in several directions: we obtain sharp regularity results (for the linear model 
problems), we can handle all the possible combinations of singular limits, and we 
obtain convergence in the best possible regularity classes. 

Our approach relies on the powerful theory of maximal L p -regularity, "H°°-func- 
tional calculus, and 7?.-boundedness, see for instance [2, 8] for a systematic intro- 
duction. 
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2. Maximal regularity 

First let us introduce suitable function spaces. Let il C R m be open and X be 
an arbitrary Banach space. By L p (il;X) and Hp(£l;X), for 1 < p < oo and 
s £ R, we denote the X-valued Lebegue and the Bessel potential space of order s, 
respectively. We will also frequently make use of the fractional Sobolev-Slobodeckij 
spaces Wp(Q; X), 1 <p < oo, s e I\Z, with norm 

\\9\\w S (n;X) = H5ll W rW ( n ! X) + X) [J Q J n \ x ^y\nMs~[s])p X dxd y ) > ( 14 ) 

where [s] denotes the largest integer smaller than s. Let T € (0, oo] and J = (0, T). 
We set 

{{u G W;{J,X) : u(0) = u'(0) = . . . = w^)(0) = 0}, 
if fc+I< s <fc + l + I, fceNU{0}, 
% S (J,X), if s<|. 

The spaces oHp(J,X) are defined analogously. Here we remind that Hh = Wp for 
k e Z and 1 < p < oo, and that W* = B° p for s e R \ Z. We refer to [14, 15] for 
more information. 

Before turning to the proofs of our main results, we add the following remarks 
on the linear two-phase Stefan problem (1) and the particularly chosen extension 
Pe determined by equation (2). 

Remark 2. (a) (l)-(2) constitutes a coupled system of equations, with the func- 
tions (v, p, Pe) to be determined. We will in the sequel often just refer to a solution 
(v,p) of (1) with the understanding that the function pe also has to be determined. 

(b) Suppose p e Wp~ 1/2p {J, L p (M. n ))nL p {J, Wp~ 1/p (M. n )) and Po E Wp~ 3/p {W l ) is 
given such that p(0) = pq. Then the diffusion equation (2) admits a unique solution 

pe e W$(J,W(R n+1 )) n L p (j,w^R n+1 )). 

This follows, for instance, from [5, Proposition 5.1], thanks to 

e -\ y \ ( l-A x )i ^ g W ^2/ P( ^n+ly 

(c) The solution ■) of equation (2) provides an extension of p(t, ■) to R ra+1 . 
We should remark that there are many possibilities to define such an extension. 
The chosen one is the most convenient for our purposes. We also remark that we 
have great freedom for the extension of po. 

Let T £ (0, oo] and set J = (0, T). By Ft wc always mean the space of given 
data (f,g,h,v ,p ), i.e., F T is given by 



¥ T = ¥ T x¥ T x¥ T x¥ T x ¥ T (5, a), 
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where 





= L p (J,L p (R n+1 )), 




F| 


= W / p _1/ ' 2p (J ! Lp(K f! 


))nL p (j,w p 2 -^(r)), 


F| 


= ^/2-l/2 P(J)ip(]| 


r))nL p (j,w^- 1/p (R n )) 




= % 2 - 2/p (IR n+1 ) 




Ff^a) 


Txr2-2/p+sg(«+<r)(2- 





Analogously, we denote by Ft the space of the solution (v,p,pe)- As was already 
pointed out in the introduction, we have, depending on the values of 5 and cr, four 
different types of spaces. For this reason we set 

E T (5,a) = E^ x E|(<5, cr) x E^ (6,a>0), 

with 

E^ = W^J,LP(W n+1 ))nL p (J,W^(M. n+1 )), 

and with E^,(<5, cr) as defined in (4) and equipped with the parameter dependent 
norm given in (5). Note that then the norm in Et(<5, cr) is given by 

\\(v,p,pe)\\e t (8,<t) = II(^P,P£;)IIe t (o,o) + £|MIe§,(i,o) + ^IHIe|.(o,i) 

for (v, p, pe) € Ft (S, a). For fixed <5, a > by interpolation it can be shown that 

E 2 (S, a) = W p ~ 1 / 2p (J, L p (R n )) n L P (J, W p i ~ 1 / p (M n )) 

in the sense of isomorphisms. We remark that Ey(c>, a) is the correct regularity class 
for the free surface if both, surface tension and kinetic undercooling are present. 
The space E^(0,ct) or Ey(<5, 0) is the proper class if just surface tension or just 
kinetic undercooling, respesctively, is present. Finally, E^(0, 0) is the correct class 
if both of them are missing, i.e., E^(0,0) is the regularity class in the situation of 
the classical Stefan problem. 

The corresponding spaces with zero time trace at the origin are denoted by oF^, 
oR}., oE 2 -((), cr). and so on, that is, 

F| = ^ p 1 - 1 / 2p (J,L p (R"))ni p (J,^ 2 - 1 /P(R")) or 

E^ = W^J,LP(R n+1 ))nL p (J,W 2 (R n+1 )), 

for instance. Moreover, we set 

F T := F^ x F| x F^, 

E T (5, cr) := Et x E|((5,cr) x E^. 

2.1. Zero time traces. We will first consider the special case that 

(h(0),g(0),v o ,p ) = (0,0,0,0). 
This allows us to derive an explicit representation for the solution of (l)-(2). 

Theorem 2.1. Let p e (3, oo), T, R > 0, < S, a < R, and set J = (0, T). Suppose 
that 

(f,9,h)€ ¥ T 

and that the function a ~ a(S,a) satisfies the conditions in (3). Then there is a 
unique solution 

(v,p,p E ) = (v»,p»,p E )£ E T (5,o-) 
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°f (!H 2 ) satisfying 

IIKp,Pb)|| e t (5,<7) < C||(/,3,/i)|| f t (15) 

with C > independent of the data, the parameters (S, a) G [0, R] 2 , and T £ (0, To] 
for fixed Tq > . 

Proof, (i) In order to be able to apply the Laplace transform in t, we consider the 
modified set of equations 

(d t + K- cA)u = f in (0, oo) x R n+1 , 
7U ± — aA x rj + 8(d t + k)t] = g on (0, oo) X R", 
(d t + K)r) + [crydyfu - arj E )\ = h on (0, oo) x R™, (16) 
u(0) = inK" +1 , 
?7(0) = hiR™, 



and 



(dt + K- cA)n E = in (0, oo) x R n+1 

vrfe = ?7 on(0,oo)xM" ! (17) 



Ve(0) = 



m 



for the unknown functions (u, 77,77s) and for a fixed number k > 1 to be chosen 
later. We claim that system (16)-(17) admits for each (f,g,h) € o^oo a unique 
solution 

(77, 77, 77s) G oEoo (<5,cr) 
satisfying inequality (15) in the corresponding norms for T = 00. 

(ii) In the following, the symbol " denotes the Laplace transform w.r.t. t combined 
with the Fourier transform w.r.t. the tangential space variable x. Applying the two 
transforms to equation (17) yields 

(u 2 - cdl)jm(y) = 0, y 6R, 

77J3 (0) = 7), 

where we set 



oj = u(\, \t\,y) = y\+K+c( y m 2 , 
w± = w±(A,iei) = x/a + c± iei 2 - 

with c(y) = c± for (±77) > 0. Equation (18) can readily be solved to the result 

77l(7/) - e"^ |y| 77. (19) 

Next, applying the transforms to (16) we obtain 

(uj 2 - cd 2 )u(y) = f(y), t/GR, 
u ± (0) + a\^\ 2 fj + S(X + K)fj = g, (20) 
(\ + K)f)+lcd y (u-arj£)(Q)] = h. 

By employing the fundamental solution 

k±(y,s) := — ( e -"±\ys\/VZ± _ e -"±fe+s)/V5±) 5 y,s>0 
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of the operator (u>± — c±dy), we make for vr 1 the ansatz 

/>oo 

u+(y) = / k + (y,s)f + (s)ds-e"+y^{<j\Z\ 2 fi + 6(\ + K)fi-g), y > 0, 
Jo 

/>oo 

u-(y) = / k-(-y,s)f-(-s)ds-e u -y^( (J \tffi + 5(\ + K)fi-g), y<0. 
Jo 

(21) 

A simple computation shows that 

d v u+(0) = — / C - u+s/ ^f + (s)ds + -J-(a\^\ 2 fj + S(X + K)fj-g) and 
c+ Jo Jc^ 



d y u-(0) = / c-"- s /^f-(- s )ds-—(a\tff ] + 5(\ + K )f 1 -g). 

c_ ./o ^/cT 



Inserting this and the fact that d y r]E (0) = ^F-^=r) in the third line of (20) yields 



m 



(22) 



with 

m(A,|€|) = A + K +(a|e| 2 +(5(A + K ))(VST W+ (A,|el) + V^^-(A,|CI)) 

+ a + V^ W+ (A,|CI)+«-V^^-(A,|ei)- (23) 

(iii) In order to show the claimed regularity for the Laplace Fourier inverse of the 
representation (u, fj) we first show regularity properties of the symbols involved. To 
this end let us introduce the operators that correspond to the time derivative and 
the Laplacian in tangential direction. Let r, s > and 

J,/Ce {H, W}. 

Then by JC s p we either mean the space HL* or the space Wp . On the space o-T 7 ^ (R+ , (1 
we define 

Gu = d t u, ue®(G)= F r p +1 (R + ,lC;(R n )), (24) 

and 

D„u = -An ue 9(D n ) = o^(K+,/C; +2 (R")), 

that is, D n denotes the canonical extension to o^p(R+,Kp(l n )) of —A in K.p(M. n ). 
Note that 

GeKF( j;(R + ,K;(l"))) with =7r/2 (25) 

and 

D„ e7eW°°(o^;(K+,X;;(R™))) with ^°°=0, (26) 

i.e. both, G and D n admit an 7?.-bounded "H°°-calculus with 'ft'H^-angle 4>q'°° = 
ir/2 and <f> D '°° = 0, respectively. Recall that an operator A admits an 7£-bounded 
"H^-calculus with 7?.% 00 -angle <^'°°, if it admits a bounded "H^-calculus and if 



K({h(A) : h G #°°(£ ), \\h\loc < 1}) < oo 

for each <p > <^'°° , where TZ(T) denotes the 7?.-bound of an operator family T C 
C(X) for a Banach space X, see [2, 8] for additional information. 
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The inverse transform of the occuring symbols can formally be regarded as func- 
tions of G and D n . We first consider the symbol u>±. The corresponding operator 
is formally given by 

F ± = (G + K + c ± D n ) 1 / 2 . (27) 
Lemma 2.2. Let 1 < p < oo and r, s > 0. Then we have that 

F± : 9{F±) -^oF r p (R+,IC s p (R n )) 

with 

@(f±) = o-f; +1/2 (k+ , fc; (R n ) ) n o r v (r + , rc; +1 (R n )) , 

is closed and invertible, where we set J- = H in case 2r G N. 

Proof. The assertion follows from [9, Proposition 2.9 and Lemma 3.1]. □ 

Next we show closedness and invcrtibility of the operator 

L := G + K + (aD n + 6(G + n)) (y/c+F + + y/cZFJ) + a +y /c^F + + a-^[cZF~ , (28) 

associated with the symbol m introduced in (23), in the space oJ r ^(M. + ,K.p(M. n )). 
We will prove invcrtibility of L and derive uniform estimates with respect to the 
parameters (8, a) in various adapted norms. In view of (25), (26), and by the 
Theorem of Kalton and Weis [7, Theorem 4.4] it essentially remains to show the 
holomorphy and the boundedness of the symbols regarded as functions of A and |£| 2 
on certain complex sectors. 

In order to obtain these estimates, the following simple lemma will be useful. 

Lemma 2.3. Let G C C" be a domain. Let fi,f% '■ G — ^ C be functions such that 
fi(z) for z G G. Then the following statements are equivalent: 



ft) -1 t t(G). 
h 

(ii) There exists a cq > such that 

\fi(z) + f2(z)\>c (\f 1 (z)\ + \f 2 (z)\), zeG. 

Proof. We set 

g.G-^-K, g{Z) — . ... . . .. , z G G, 

l/lW | + |/2(2:)| 

which is a well defined function. Observe that (ii) is equivalent to saying that 
^ g{G). By contradiction arguments it is not difficult to show that this relation 
is equivalent to condition (i). □ 



Remark 3. The assumption fi(z) ^ for z 6 G is just for technical reasons and 
can be removed. 

Now we prove closedness and invcrtibility of L. 

Proposition 1. Let 1 < p < 00, r, s > 0, R > 0, (8, a) G [0,i?] 2 , and F,JC G 
{H, W}. Suppose that a is a function satisfying condition (3). Then there is a 
number n > I such that 

9{L) = ^ +1+sg(5)/2 (R+,^(R n ))n ^p +1 (lR+,/c; +sg(5) (R™)) 

n o^I +1/2 (M + , /c; +2sg(ff) (K n )) n oK(^+, /q; +1+2ss(ff) (»")) 



10 



JAN PRuSS, JuRGEN SAAL, AND GIERI SIMONETT 



and L : 5${L) — > n J :7 p (R + , K.p(R n )) is invertible. Furthermore, 

aWD^G+lf^L-^ + aWD^L- 1 ^ 

+ 5\\(G + lf' 2 L~% + SWDl/^G + l)L-% + \\L-% < C 

with C > independent of (5, a) G [0,R] 2 , where |j • || denotes the norm in 

^f(o^(R + ,^(R"))), 

and || • ||i the norm in 

( t; (k+ , fc; (R" ) ) , Q r p +1 (R+ , fc; (k- ) ) n r r p (r+ , k°+ 1 (R" ) )) . 

Proof. Let tpo G (0,7r/2) and tp G (0, tpo). By a compactness and homogeneity 
argument it easily follows that 



|w±(A,2)| = IV A + K + c ± z \ 

> c [^\X\ + ^+c±VW\) (29) 

for all (A, z, re) G T, 7r ^ Va x S v x [1, oo) and some cq > 0. 

In the following we let tpo G (n/3,n/2) and tp G (0, tpo — 7r/3). Note that by 
condition (3) on a there exist 5* , <j* > and M, Co > such that 

a±(S,<r)>co {(6, a) G [0,6*] x [0,<r*]) (30) 

and 

|a ± (<5,a)| <M ((<S, a) G [0, i?] x [0, i?]). (31) 

First assume that (30) is satisfied, i.e., that (6, a) G [0,5*] x [0,<r*]. Let m be as 
given in (23). We consider the function 

/ : S w _ yo x x [0, 5*] x [0, a*} x [1, oo) -> C, 

(A,z,5, cr, re) h-> f(\,z,5,o-, re) := m(A, z) := /i(A, z, <5, cr, re) + / 2 (A,z,<5,cr,re), 
with 

fx{\,z,a,8,n) := (A + re) [<5( v / c^'aj + (A, z) + y^cTu^A, z)) + l] , 
f2(X,z,a,6,K) := m(A, z) — /i (A, z, cr, 5, re) 

= ( v ^jTw + (A,z) + ^/cTw_(A,z)) 

+ a + (<5, cr) v /c^'aj + (A, z) + a~(S, fr) v / cTw-(A, z). 

Note that by our choice of the angle tp for (A, z, <5, cr, re) G S w _ vo x E v x [0, J*] x 
[0,cr*] x [l,oo) with argA > there exists an e > such that 

7T — 9?0 /T 3<£> 3</3o I" 

- po > — 2 h ^ - ar S^V A + re + c ± z > — — > — + - + e, 

if cr > 0, and that 

K — PO /T f 

— ^ — > ar S V A + « + c±z > --. 

By these two estimates we see that in any case we obtain 

3(?r - go) , . p 
^ > ar 8/i(A,z,0,CT,re) > --. 

and 

7T - <^o > /a (A, Z, d, cr, re) > — — > — + - + e. 
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Consequently, 

2tt w f 2 (X,z,5,a,n) 3ip ir 3(tt - tp ) 

— >n-ip + ->axg — > — 7r+o +£ o 

3 2 j i (A, z, o, cr, k) 2 2 2 

A similar argument holds for the case that (A, z, S, a, k) G S t _ vo x S v x [0,(5*] 
[0,cr*] x [l,oo) with argA < 0. Here we obtain 

2tt / 2 (A,2,(5,<7,k) 

— 5- < arg — — < 7T - e. 

3 /i(A, 2,d, cr, K) 

This implies that 



-1 £ S 7r _ E D ^ x E v X [0,5*] X [0,(7*] X [l,oo)) 

Jl 

Lemma 2.3 now yields the existence of a c\ > such that 

|/i(A, z, S, cr, k) + f 2 (X,z, (5,cr, k)\ > ci(|/i(A,2;,(5,(7, k)| + \f 2 (X, z,5,<j, k)\) 

for all (A, z, S, a, n) G E„._ vo x x [0, (5*] x [0, a*] x [1, oo). An iterative application 
of Lemma 2.3 on the summands of f\ and f 2 and an application of inequality (29) 
then result in 

\f(X,z,a, 5,k)\ 

> c 2 ||A| + K + a\z\ L/\X\ + y/K+ y/c+\z\ + V c -\ z \) 
+ 8(\X\ + k) + + v^TN + v^R) 
+ a + (y\X\ + \/k + Vc+R) + a - (V\M + Vk+ v/c-I^Q } , 

for all (A, z, S, a, k) G Ett- Vo x T, v x [0,6*] x [0,cr*] x [l,oo). This implies that the 
functions 



1 A + K 

mo ■= j, mi := - , 



v/i azy/X + k 

m 2 : - m 3 := — 



az 3 ' 2 5(X + k) 3 / 2 

mi ■= — — , m 5 := ■ 



I 



m e := 



6(X + k)^/z 



are uniformly bounded on S w _ vo x E v x [0,(5*] x [0,cr*] x [l,oo). 
Now consider the cases i?><5><5*>0or_R>cr>(7*>0. We set 

g(X,z,5,a,n) := f(X,z,S,a, k) - a+(S, a)y/c^uj+(X, z) - a-(5,<r)y/cIw-(\, z). 

The argumentation above shows that 

1 X + n azy/XT^ az 3 / 2 S(X + k) 3 / 2 S(X + K)y/z 

7 7 7 7 7 

5 5 9 9 9 9 

are still uniformly bounded functions and this even on Yj 7T ^ Vo x E v x [0, R] 2 x [1, oo). 
The aim now is to show that the term a+(6, a)^J~c+uj+(X, z) + a~(8, a)^fc~uj—{X, z) 
can be regarded as a perturbation of g, if k is assumed to be large enough. Indeed, 
if 5 > <5* > 0, by using (31) we can estimate 



a±(S,a)JcIu± 



< 



CM 



5*\X + k\ 
C 



S(X + k)u>± 



C 
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for (A, z, S, a, k) e E n -tp x £ v x [<J*,jR] x [0,R] x [1, oo). On the other hand, if 
cr > a* > 0, we deduce by virtue of (29) that 



ft 



< 



< 



< 



CM 

c 

C_ 



A + K + c±z 



A + K 



a*\\f\ + k\ 



for (\,z,5,ct,k) e E T _ V0 x S v x [0,R] x [a* ,R] x [l,oo). Hence, for fixed n chosen 
large enough we see that we can achieve 

a+(6, a)y/c^uj+ + a-(5, a)y/clLJ- . 1 



< 



to be valid for (A, z, 5, a) € Sjt-^o x x x [0, -R] or (A, z, 6, a) £ 

E v x [0, -R] x [cr*,i?]. Thus, we may represent 1/f as 

1 1 / a+(<5, a)y/c+(*)+ + a— (5, a)y/cluj- 
1 9 V 3 
and therefore the functions mo, . . . , m6 are uniformly bounded for all (A, z, S, a) £ 

^'7T ~pa X ^"V X "0* R" • 

The remaining argumentation is now analogous to the proof of Lemma 2.2. Em- 
ploying (26) we obtain 

Ai,<5,o-)||jSf( ^(M + x;(K"))) : (A, 8, er) e x [0,i?] 2 |^j 

for j = 0,1, ... ,6. Consequently, 

||m J (G, J D„,<5,a)|| w;(R+ ^ (Rre))) < C ((S,a) e [0,i?] 2 ), 

by virtue of (25) and [7, Theorem 4.4]. The invertibility of the operators 

(G+l) 1 / 2 : ^ +1 / 2 (R + ,/C;(R"))^o^(R+,/C;(M")), 

DU 2 + l : ^(K + ,/C; +1 (K"))^o^(K + ,^(R")), 

(see for instance Proposition 2.9 and Lemma 3.1 in [9]) then yields the assertion, 
since L^ 1 = mo(G, D n , a, 5), and by employing the fact that h i— > /i(G) is an algebra 
homomorphism from H°°{Y,^, PQ ,X G {X)) into Jgf(AT) for X = -7>(R+, £p(R n )) 
and where 

JT G (X) := {B e if (X) : B(ji - G)" 1 = (/i - G^B, p g p(G)}. 

□ 

(iv) We turn to the proof of the corresponding regularity assertions in Theorem 2.1 
for (u,T],rjE)- According to the results in [5, pages 15-16], 



e -F + s/ V c Tf+{s)ds 6 ^ 

By the same arguments we also have 

f°c- F - s /^r(- s )d,se FL 



/+el# +) l p (lf)). (32) 



re£ P (M+,I P (M! +1 )). (33) 



SINGULAR LIMITS FOR THE TWO-PHASE STEFAN 



13 



Next, note that by Lemma 2.2 we have that 

F ± eIsom( F^, F^). (34) 

Indeed, we obtain 

F^Wl) = ^ 1/2p (K + ,L p (]R n ))no^y 2 - 1 / 2p (]R + ,% 1 (R")) 

n oH^/ 2 (M + ,Wp~ 1 ^ p (R n )) n L p (m + ,w^ p (R n )) 



? 2 

■ OO 1 



by virtue of the embedding 

oF 2 ^ ^ 2 - 1 / 2p (M + ,^ p 1 (R' 1 ))no^ /2 (M+,% 1 - 1/p (M' 1 )), 

which is a consequence of the mixed derivative theorem. Thus all the terms inside 
the brackets on the right hand side of (22) belong to the space o^L- I n t ne same 
way as we clarified the invertibility of F± : oF 2 ^ — > oF^, by applying Lemma 2.2, we 
can see that L : oE 2 c ((5, a) —> q¥^ >0 is invertible by an application of Proposition 1. 
For instance, if 6, a > 0, this follows from the embedding 

oE^a) ^ o^/ 2 (R+, W^(R n )) n ^- 1/2,, (i +) ^(l")), 

which is again a consequence of the mixed derivative theorem. Furthermore, Propo- 
sition 1 implies the estimate 

l|i _1 !l^(oF^,oE^(0,0)) + <5||i _1 ||^f(oF3 c ,oE2 (l,0)) + &\\L~ II^CoF^.oK^CO.l)) < C 

for < 6, a < R. Altogether this gives us 

IMU^.a) ^ C (ll/lk + HflHoFL + IHUfJ (35) 

for ((5, a) e [0,i?] 2 , which yields the desired regularity for rj. Observe that u now 
can be regarded as the solution of the diffusion equation 

(d t + K- cA)u = f in (0, oo) x W l+1 

7u ± = g + aA x rj — S(dt + K)r] on (0,oo) x R n , 
u(0) = inR n+1 . 

A trivial but important observation now is that this equation itself does not depend 
on S and c, but only the data. Therefore also the corresponding solution operator 
is independent of 5 and a. By well-known results (see e.g. [5, Proposition 5.1]) and 
in view of (35) we obtain 

IMIoE^ < C (II/IIf^ + llslloF^ + (S hlloE^(l,0) + ^hll E^(0,l)) 

< C\\(J,g,h)\\ oTae (0<S,a<R). 
Similarly we can proceed for r\E- Since it satisfies equation (17), we deduce 

IMU^ < ClhlloF^. 
By virtue of E^(0,0) '—t oF 2 ^ and again (35) we conclude that 

II^HoEJ. <C||(/,5^)lloFo. (0<<y,ff< Ji). 

(v) Let T > be fixed, and let J := (0, T) with T < T . We set 
TZj: Ft -> oFoo, 
(f,g,h) i y (e- Kt (£jf),e- Kt (£jg),e- Kt (Sjh)), 
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where £j is defined as 

( u{t) if < t < T, 

£ju(t) := £j. r u(t) := I u(2T -t) if T < t < 2T, 

[ if 2T < t. 

It follows from [10, Proposition 6.1] and the fact 

\\(e- Kt (£jf),(e- Kt (Sjg),e- Kt (£jh))\\ oroo < \\e- Kt \\ BVCl(K+) \\(£jf,£jg,£jh)\\ oVoo 

that there exists a positive constant cq = cq(Tq) such that 

m(f,9,h)\\ o¥aa <c \\(f,g,h)\\ oWT ((f,g,h)e W T ) (37) 

for any interval J = (0, T) with T < To- 
Let (u,r],r]E) & oEoo(i5, c) be the solution of (16)-(17), with (f,g,h) replaced by 
{7lj(f, g,h)), whose existence has been established in steps (i)-(iv) of the proof. We 
note that 

II {u,V, f/B)|| E oo (5,<T) < K\\1Zj(f,g,h)\\ oVco 
< Kco\\(f,9,h)\\ w T 

for any (f,g,h) 6 Ft, < 5, a < R, and any interval J = (0,T) with T < T , 
where K is a universal constant. Now, let 

{v,p,p E ) := (R J {e Kt u),'R.j{e Kt i 1 ),'Rj{e Kt i lE )) 

where IZj denotes the restriction operator, defined by IZjw := w\j for w : R+ — > X. 
Then it is easy to verify that 

(v,p,p E ) £ E T (6,a), (v,p,p E ) solves (l)-(2) (38) 

and that there is a constant M = M(T ) such that 

\\(v,P,PE)\\ E T (6,a) < M||(/,g,/l)|| oFT 

for < S, a < R, and T < Tq. Finally, uniqueness follows by a direct calculation 
which is straight forward and therefore omitted here. This completes the proof. □ 

We proceed with convergence results for the case of zero time traces. To indicate 
the dependence on the parameters 6 and a we label from now on the corresponding 
functions and operators by /i, as e.g. L^, v^, where p = (5, a). 

Corollary 1. Let 1 < p < oo, R > 0, < S < S < R, and < a Q < a < R. 

Suppose that a is a function satisfying the conditions in (3), and let be the 
operator defined in (28) corresponding to the parameter p := (5, a). Then we have 

(5-5 )L- 1 ^0 stronlyin if ( F^, E^(1, 0)), (39) 
(a - (Jo)^ 1 -> stronlyin ££ ( F^, E^(0, 1)), (40) 

and 

L-^L-l stronlyin -§f( F^, oC(Mo)), (41) 
as p — >• po, where po = (<5q, ctq). 
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Proof. As pointed out in part (iv) of the proof of Theorem 2.1 the domain of the 
operator F+ in qF^ is oF^. This implies that 

9{F%) ^ 1/2p (M+, L P (W 1 )) n L P (R + , W^(R n )) = E^(1, 1). 

Now pick / <E @(F+). From Proposition 1 we infer that 

II^IUcoFL.oiP.Co.o)) < C (/i G [0,i?] 2 ). (42) 

This yields 

||(5 - S )L^f\\ oEUm < C(S S ) (|| (G + «) 3/2 i; 1 /ll w v- 1 / 2 , (H+>Lp(R „ )) 

+ ll( G + K ) i ^ 1 /ll Lp(R+!V l / 2-l/P (R „ )) ) 

< c(<y-«o)||/|| E 3o (i, 1 ) 

-> (n^no). 

Since @(F+) is dense in oF^, the uniform boundedness of H-ku^ljSfdjFS , e 2 m 
M € [0, R] 2 (which yields uniform boundedness of (S — ^o)!!-^ 1 !! if(oF^ , E^ (i,o)) f° r 
5 £ [So,R], a £ [0, R)) implies (39). In a very similar way (40) can be proved. In 
order to see (41) we write 

L »o ~ L v 1 = L vl ( L n ~ L »o ) L~ 1 

= {((5 - 5 a )(G + *) + (*- *o)A0 {^F + + ^fc~F_) } L^ 1 

+ ^ {(<*+(**) - a+(Mo))^F+ + (a-(fi) - a^))^IF^}L^\ 

In view of £ (oF^, oE 2 (mo)) this representation shows that (41) is obtained 
as a consequence of (39)- (40), and (42) in conjunction with the continuity of a± . □ 

Based on this result we will now prove convergence of solutions of problem (l)-(2). 

Theorem 2.4. Let 3 < p < oo, R, T > 0, < 5 < 5 < R, and < a < a < R. 
Suppose that a is a function satisfying the conditions in (3) and that 

((r,g^h^)^ [So , R]x[a0iB] c F T . 

Furthermore, denote by {v^, p^) the unique solution of (l)-(2) whose existence 
is established in Theorem 2. 1 and that corresponds to the parameter p = (5, a) . 
Then, if 

(r,g^h^^(f^,g^,h^) m F T (m ^ Mo), (43) 

we have that 

(v*,pP,p%)^W«,pi«,p%) tn E T ( Mo ) (m^Mo), (44) 

where fio — ((5o,oo)- In particular, if 

S- 1 :{f,g,h)^(v^p\p^) 

denotes the solution operator to system (1), we have that 

S ^ 1 -> S ^o strongly in if ( F X , E^ X E|(mo) X E^) (fj, -> Mo)- (45) 

Proof. In view of the arguments in part (v) of the proof of Theorem 2.1 the solution 
(w M , p M , p^) can be represented by 

(v",p",p£) := (T^V^^V),^^)), (46) 

where 7£j denotes the restriction operator and (m m , rf , rfe) is the solution of (16)- 
(17) with right hand side (7?.j(/ M , g^, h 11 )) and IZj as defined in (36). Hence we see 
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that it suffices to prove convergence for the vector (it^, rf , 77^.). Clear iy, (43) implies 
that 

(Tlj(r,g^h^)^(n^r,g^,h^)) in Foo (ji -> Mo). 

Therefore, and for simplicity, wc simlpy write {f^,g^,h^) for the data instead of 
(7£j(/ M , g^jh^)) in the remaining part of the proof. 
Next, recall from (22) that t/ 1 is given by 

rf=L- x l» (/xG [* ,oo) x [0-0,00)) 

with 

/* OO />OG 

r=^-/ e- F + s ^{f) + (s)ds- c- F - s /v^(/^)-(-s)ds 
Jo Jo 

According to (34) we know that F± £ Isom(oF^ o; oF^,). This fact and relations (32) 
and (33) then imply, by virtue of assumption (43), that 

r - '"iioFL < c (wr - n^ T + - <? mo iui + w - ^°iu^) 

-> (n^no). 

By the uniform boundcdness of H-k^ 1 !! ^( w 3 , e 2 O)) in /j G [$0;-R] x [<7o,-R] (see 
Proposition 1) and because (S — So)^ 1 —> strongly in Jzf (oF^, oE^(l, 0)) and 
(<t - aajL' 1 — > strongly in Jz?( F^, oEt(Oi 1)) ( see Corollary 1) this results in 

(8-So)if^0 in E^(1,0) (47) 

and 

(a-<7o)<^0 in E»(0,1). (48) 

Now, denote by 

U» : K,^)^ (f^g^h^ 

the operator that maps the solution to the data corresponding to system (16). From 
part (iv) of the proof of Theorem 2.1 we infer that 

U„ £ Isom( E^ x E^ (p) , qF^ ) (ji £ [S ,R]x[a ,R]). (49) 

Furthermore, observe that we have 

= U-\f»,g*,h») - U- ( ?(r° 7 g^,h^) 

- g^> + a a )A xV » -(5- 6 )(d t + n)rf 

Relation (49) applied for n = no then yields 

\\(u^,^) - (u^, V ^)\\ oKloXoMUflo) 

<cf||(^^,^)-(/*,^,ft^)|| F ao + (*-<5o)||i^|| oESo (i,o) 
+ (ff - ^o)||^|| e^(o,i) + HfJ-) - aMWlVsWo^ 
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From Theorem 2.1 we know that ||t7jj|| el, is uniformly bounded in p € Iq. Thus, 
by (3), (47), (48), and assumption (43) we conclude that 

(«",*?") ->• (u" ,*?" ) in oE^XoE^Oio) (/" ->• Mo)- 

The convergence of 7?^. is easily obtained as a consequence of the convergence of 
if-. Recall that rf^ is the solution of (2) with p replaced by rf . Denote by T the 
solution operator of this diffusion equation which is obviously independent of p. 
Then by [5, Proposition 5.1] we obtain 

K-P£°IU L = HT(0,^-rr,0)|U L 

< CIK-^llo^Co.o) 

-> (^->Mo), (50) 

by the just established convergence of if . Representation (46) then implies (44). 
Obviously (44) is still true for fixed data, i.e., if 

(F,g^hn = (f,g,h) e Ft (p G [5o,R] x [o- ,iJ]). 
Hence (45) readily follows from (44). □ 

2.2. Inhomogeneous time traces. Next we consider the fully inhomogeneous 
system (1)— (2) and we will prove Theorem 1.1. By introducing appropriate auxiliary 
functions, we will reduce this problem to the situation of Theorem 2.1. 



Proof, (of Theorem 1.1.) If S = a = this result is proved in [10, Theorem 3.4] 1 . 
So, we may assume that 6 > or a > which implies that po G Wp 3 ^ p (R n ). 
Furthermore, it follows from the trace results in [3] that the conditions listed in 
(6)-(8) are necessary. 

Suppose we had a solution (v,p,pe) of (l)-(2) as claimed in the statement of 
Theorem 1.1. Let v\ be the solution of the two-phase diffusion equation 

(dt-cA)ui = / in JxE n+1 , 

jvf = g + e'^-^C onJxr, (51) 
ui(0) = w inM" +1 , 
with 

C:= 7*0-5(0). (52) 
Observe that by compatibility assumption (7) we have 

( = (aA x p-6d t p)\ t = . (53) 

Next let pi be an extension function so that 

( Pl (0),d tPl (0)) := (p 0) /i(0) - [ndyivo-ae-^-^po)}) , (54) 

as constructed in Lemma 3.2, and let p\^E be the solution of (2), with p replaced by 
p\. For the solvability of (51) and the existence of p\ we have to check the required 
regularity and compatibility conditions for the data. By construction we have that 
<?(0) + £ = 7«o and by the regularity assumptions on g and vq we deduce 

C = jvo - g(0) G W^ 3/P (^ n ), 



i Actually with g = 0. But by obvious changes in the proof one can obtain the result also for 

o + g e f£ . 
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hence that 

e -(l-A.)t£ g F 2 _ (55) 

Then it follows from [5, Proposition 5.1] that there is a unique solution Vi € of 
(51). Furthermore, if 6 > 0, we may use compatibility condition (7) to obtain that 

h(0) - lndy(v - ae-W 1 -**) 1 " pa)] = 1( 5 (0) - 7 i'o + aA xPo ) e % 2 - 3 ^(M"). 
If S — 0, we may impose a > which gives 

c 7 d H ae-l^ 1 - A ^ 1/2 ( oo = Tca(l - A x )^ 2 p G VK p 3 - 3/p (M") VK 2 ~ 6/p (R") 
in view of pa £ Wp 3 ^ p (W l ). Assumption (8) then implies that 

MO) - lcjd y (v - ae-MC 1 -*-) 1 " ' p )j e H/^ 6/p (l"). 

Thus, in any case we can satisfy the assumptions of Lemma 3.2 which yields the 
existence of p\ € ¥^(5, a) as claimed, and of p\^ E G E^ by virtue of Remark 2(b). 
Now we set 

(v2,P2,P2,e) = {v,p,p E ) - (v 1 ,p 1 ,p 1 , E ). 

It is clear that p2,E is the extension of p2 given by (2) with p replaced by p2- Thus, 
(«2j P2, P2,e) satisfies 

{d t -cA)v 2 = inJxi n+1 , 
fvf - (tA xP2 + 6d t p 2 = aA xPl -Sd t pi-e- ( - 1 - A ^ t C onJxi", 
dtp2 + lcyd v (v 2 - ap 2 ,E)] = h - d t pi ~ \c~jd v {v x - ap hE )j on J x W l , 
v 2 (0) = inR" +1 , 
p 2 {0) = inM n , 

(56) 

and 

(dt-cA)p 2 ,E = mJxR n+1 , 

7P± B = p 2 onJxR™, (57) 
P2,e(0) = ini" +1 . 
By construction, p 1 6 E^(<5, a), and by (55) one may readily check that 

oA xPx - Sd tPl - e -(!- A *)*C g F| 

and that 

h - d t pi ~ [cydyiyi - ap hE )j G F|- 

Thus, by Theorem 2.1 the reduced system (56)-(57) is uniquely solvable. This 
allows us to reverse the argument. In fact, since the solution v\ of (51) and the 
extension p\ depend on the data only, the right hand side of (56)-(57) so does as 
well. Theorem 2.1 now yields a unique solution (v2, P2, P2,e) G oEt(<5, c) and 

(v,P,Pe) ■= (v2,P2,P2,e) + (vi,Pi,Pi,e) (58) 

then solves the original system (l)-(2) in the reguarity classes required. It remains 
to verify estimate (9). Observe that by Theorem 2.1 we know that 

\\(V2, P2, P2,e)\\ E t (S,<j) 

< C (\\aA x p x - Sd tPl - e-^-^'ClloF* +\\h- 8 tPl - [c 7 9,(«i - a Pl , E )]\\ oK ) 
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with C > independent of S, a. By \a(p)\ < C for /i £ [0, R} 2 and the facts pointed 
out above we can continue this calculation to the result 

\\{v 2 ,P2,P2,e)\\ 

< C (a\\/\ xPl \\ ¥ 2 T + 6\\d t pi\\ F 2 T + \\(\\ w 2- 3/P{Rn) + \\h\\ ¥ 3 T 

+ \\dtpi\\v 3 T + IK - api,E\\uC) 

< G (\\(vi,Pi, Pi,e)\\e t (8m) + ||(0,3,Mo,0)||]ev(o,o)) ■ ( 59 ) 

Hence we see that it remains to derive suitable estimates for (i)i,pi, Pi,e)- Observe 
that equation (51) does not depend on 5, a. By [5, Proposition 5.1] we deduce 

IKIU < c(\\f\\ ¥lT + \\g + e-^ A ^(H + \\vo\\ r i T ) 

< C (II/IIfi, + HffllF- + II«oI!f* ) (0<5,a<R). (60) 
By the same argument we also have 

< <?(|MIe|(o,o) + \\Po\\ w *-VP (Rn) ) (0<S,a<R), (61) 

where we used Remark 2(b) and the embeddings Wp 2_2/p (R' 1 ) ^ Wp 2_3/p (M") and 
Ey(0, 0) ^ F|,. Lemma 3.2 implies for pi, 

l|Pl|| o E|(0,0) 

< C (HpoII^-Vp^) + ||^(0) - lcjd y (v ~ ae-^-^pojW^-^^) 

< c + \\h\\^ T + \\v q \\ ¥ 4 t + ||(1 - A x ) 1/2 / 9o]j| Wp i-3/P (R „ ) ) 

and 

°1|Pl|| E§,(0,l) 

< C (vWpoW^-V^ + a\\h(0) - indyivo - ae^'^poM^-e^^) 

< C {^WpoW^-Vp^) + cr||/i,(0) - [crf9yVoj\\ w 2-6/ P{Mn) J (63) 
as well as 

<5 l|Pl|| E|,(l,0) 

< C (s\\p \\ wi - 3/P{Mn) + 6\\h(0) - Icjdyivo - ae-^-^poM^^^) 

< C ( s \\Po\\ w i-i/P {Rn) + \\g\\w 2 T + INIf^, + cr ll/ ? o|| H -4-. V p (R „ ) ) , (64) 

for < 6, a < R, where we used in (64) once again compatibility condition (7). 
Inserting (62) into (61) we obtain by (60)-(64) that 



IK^i, Pi,Pi,e)\\e t (8,<7) < C (\\(f,g,h,v ,p )\\ VT(Qfi) + (8 + o-)||poll w -4-3/P ( 

+ o-|IM°) - l c ldyV ]\\ w 2- B/P(Mn) j (65) 
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for < S, a < R. Inserting (65) into (59) we can derive exactly the same estimate 
for («2, P2, P2,e)- Combining the estimates for (vi, pi, Pi,e) and (i>2, P2, P2,e) we 
finally arrive at (9) and the proof is complete. □ 

Next we prove convergence for the solutions of problem (l)-(2), that is, Theo- 
rem 1.2. 



Proof, (of Theorem 1.2.) 

We employ the decomposition 

P» = P?+P2 l 

as given in (58). We have to show that 

(i) K.pW.J^^.rf ,^) in E r ( Mo ), 

(ii) K,^,^, £ )^K°,^°,P^) in ErOxo). 

(i) We start with proving convergence of p p . This function is according to (54) an 
extension of the traces 

where we set 



(pf(0),ftpf(0)) := (p&,<#) 



(66) 



Since the extension operator in Lemma 3.2 is linear and independent of p we can 
estimate for all ji € Jq, 



where 



Hpi -pHk(Mo) < c(llp < -Po to lk(Mo) + lko-?o <0 lk( W )) 



Ff^o) := W^- 3/p (M. n ) n w^ aoK2 - 6M (m. n ) n w* s{So){2 - 3/p) (W l ). 



(67) 



It is clear by (10) that the first term on the right hand side of (67) tends to zero. 
In order to see the convergence of the second term we distinguish the three cases 
So = &o = 0, and <5o > 0,o~o > 0, and <5o = 0, cto > 0. 

The case 5a = ctq = 0: Here we have Fy(/zo) = Wp 3 ^ p (W l ) and we obtain by a 
direct estimate and (10) that 

MM p-o 1 1 

n?o - % Hw,;- 3/ "(R") 

< c (\\h» - h">\\ n + |K <°IIf* +U- pi°\\ w ^ (m ) 

-^0 (p^p a ). 

The case 6 > 0, cr > 0: Then F|(/x ) = Wp~ 3/p (R n ). In this case we can employ 
compatibility condition (7) in Theorem 1.1 which results in 



2-3/p, 



i GA0) - 7 < + aA x p%) - 1 (g»°(0) - 7 <° + o-A x p%°) 

Oo 



< C 



1 



1 



1 



1 



■Po 



A ^° 
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In view of So > observe that p£ -> p^° in F^(/i ) = Wp~ 3/p (R n ) by (10). This 
yields 

s a II P Moil i ^ °"0 \ ., W) | 

S •jHPo - Po IIfIGko) + I ^ - ^ I IlPo IIf|,(mo)- 



Ff,( Mo ) 

-> (n^fio). 

In the same way we see that the first and the second term on the right hand side 
of (68) vanish for p, — > p,Q. 

The case Sq — 0, cto > 0: Since <$ — > 0, here we cannot apply compatibility condition 
(7). This leads to condition (11) in the statement of the theorem. In fact, here we 
obtain 

H - Cllv^( K ») < c (ll^(o) - [cydX] - fc"°(o) + MtClll^-/,^) 

■ II j-L /./.on \ 

+ IlPo -Po Hm/ p 3 - 3/p (R")J ■ 

It is clear that for Co > condition (11) implies that the first term on the right 
hand side vanishes, whereas the second term tends to zero again by (10). 

Also here the convergence of p^ E follows by the convergence of pi in view of the 
fact that p^ E is the solution of (2) with p replaced by p%. If T denotes again the 
solution operator of this diffusion equation, by [5, Proposition 5.1] we obtain 

= lir(o,^-pr,e-' y ' (1 - A ^" /2 (^-^°))|| E i 

< C (JK - pi? || F| + lle-l^ 1 -**)- 172 ^ - ^)|| F 4 ) 

< c (up? - pfik(o,o) + M - 

-> (M^Mo), (69) 

by the just proved convergence of p± and (10). 

Observe that v% is, according to (51), the solution of the same diffusion equation 
with right hand side (/ M , + e _ ( 1_Aa, )*(7i>o — g^(0)), Vq) for p € I - Moreover, we 
have that 

!|e- (1 - A * )4 (7< -5 M (0) -7<° +^°(0))|| F | 
<q| 7 <-5 A1 (0)-7<°+.9 AIO (0)|| wr 3 /P(R „ ) 

Hence we obtain 

ii< - <°ii E - < c (\\f r°ii F L + ik - <°ik + - .9 wi Hf|) 

-> (/x -> /xo) 

by (10), and (i) is proved. 

(ii) Note that (v%,p2, P^e) * s * ne solution of (56)-(57). According to Theorem 2.4 
it therefore suffices to prove convergence for the corresponding data. To be precise, 
it remains to show that 

in oIt (m~>Mo), (70) 

where 



«T = aA x p1 - 5d tP » - e-d-^V, 
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and that 

fr-th" in olr (p->f*o), (71) 

where 

h" = h» - d t fi - [crydyW - a(fi) P 1 tE )j. 

First we estimate 

< c (||^ - ^°|| F s + 1|< - <°|| E ^ 

+ IKaOPi.jj - a(Mo)Pi,° B ||Ei + IK - Pi°IIe|(o,o)) , 
and we see that (71) follows from (i), (3), and (10). For g M we have 

\\r -£HIf| < C (\\5pt - 5op?\\ %W) + \\<rp>£ - ctop1 \\eI(o,i) 

+ 11^-^11^-3/^)- (72) 

By employing the convergence assumptions also here we will prove that each single 
term on the right hand side of (72) tends to zero for p — > po- In view of (52) and 
(10) the convergence of the third term in (72) is clear. The first two terms are 
more involved. In fact, this is the point where assumption (12) enters. In analogy 
to (i) wc again distinguish the three cases Sq = ctq = 0, and Sq > 0,<7o > 0, and 
S = 0, a Q > 0. 

The case So = = 0: Note that in this case condition (7) for po turns into 

7ttf°-<r(0)=0. (73) 

By using this fact, Lemma 3.2, (7) for p, and recalling that q$ still denotes the 
function defined in (66) we obtain 

II^PiIIe§,(i,o) 

< C (5\\p»\\ w} - 3/V{Mn) + ||<f (0) - 7 < - <T>(0) + 7 <° + aA x p^\\ w ,- 3/P{Mn) ) 

<c[(5 + <j)M\\ wt ^ {Rn) + \\g" - ^|| F = + IK - <°lk) ■ 

In view of (10) and (12) we conclude that 

II<KIIe§,(i,o) -> (m~>Mo)- 
For the second term in (72) Lemma 3.2 yields 

lkPillE=(o,i) < C (^llPollw^-s/p^n) + f7 lkoll Wp 2 - 6 ^(R«)) 

< C ( (T ll^o!lw p 4 - 3 ^(R") +°1l /lAi (°) - l c "/ d y v ol\\w*- 6 'v (r»)) • 
Hence, if 5 > 0, it follows 

\Wpi ||e§.(o,i) (/«->• Mo) (74) 

by (11) and (12). If 6 = 0, we have erA^ = 7«o ~ M ( )- This y ields 

IkPollv^-^CR") - C (^ll^ollvi^-^CR") + H crA ^/ , oll W ' p 2 - 3/p (R")) 

< ^^ll^llwr^^ + ll^-^lk + ll^-^lk) , 
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where we used again (73). Observe that \\pfi || w i-2/ P( ^ n ^ is uniformly bounded in 
p € /o by assumption (10). Thus, in this case (74) is obtained as a consequence of 
(7), (10), and (11). 

The case Sq = 0, cto > 0: Here we have 

-yvp-aAxfC =^°(0). (75) 
In a similar way as in the previous case we deduce, if 5 > 0, that 

II<MlIIe?,(i,o) <C UWPqWw^-Vp^ + IkPo - <roPo°\\ w 4-3/p ( j in) 

x (76) 

+ II^-^ IIf|. + II<-< II^)- 

Note that in the case ao > we also have that 

fi^fC in ¥ 5 T ( Mo ) = <- 3/p (K") (M^Mo). 

By this fact it is easy to see that the first two terms in (76) vanish for (p — > /io), 
whereas the convergence of the last two terms follows again by (10). That the 
second term in (72) tends to zero here follows easily from the inequality 

hPi -0o/?i o ||e=(o,i) < — IK -Pi°IIe=(o,<t ) + - — —\\Pi° IIeS,(0,o- ) 

CO °0 

and the convergence of pt^ in E§,(/i ) = E§,(0, oo) proved in (i). Observe that the 
last argument also implies convergence for the case 8 = 0, since then the first term 
in (72) vanishes completely. 

The case So > 0,ao > 0: Here the convergence of the first term in (72) follows 
completely analogous to the convergence of the second term in the previous case. If 
we suppose that also <7o > the convergence of the second term in (72) follows by 
the same argument. In the case that oq = also here an application of Lemma 3.2 
implies 

IkPi 11^,(0,1) ^ c(^IIPollw p J - 3 /"(R«) +cr ll9olli Vp 2 - 6 ^(R")) 

Moreover, we still have Wj,{po) = Wp~ 3 ^ p (R n ), which implies Po II vk 4 " 3/:p (K" ) — * ^ 
for p — > p by (10). For the second term on the right hand side of the above 
inequality note that from the case So > 0, a > in (i) we know that <7q — ¥ q^ in 
W / p _3 ^ p (lR n ). This implies that this term vanishes as well for p — > /io- Hence also 
in this case we have that 

HlPi He*(o,i) -> (p^po). 
The three cases together show that 

(0,ff", &")-»■ (O,^,^) in F T (M^Mo), 
and therefore Theorem 2.4 implies (ii). □ 

3. Appendix 

The reduction of problem (l)-(2) to the case of vanishing traces in the proof of 
Theorem 1.1 was based on the following two results. Observe that the assertions in 
Lemma 3.1 follow directly from the general trace result [4, Theorem 4.5]. However, 
for the sake of completeness and for a better understanding of the proof of subse- 
quent Lemma 3.2 we give its proof here. In the following we adopt the notation of 
Section 2.2. 
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Lemma 3.1. Let 1 < p < oo, T e (0,oo], and J = (0,T). 

(i) For each t]q G Wp 3 ^ p (]R n ) there exists an extension 

such that 771 (0) = <7o an d, if P > 3, also that 9 t 771(0) = 0. 
(ii) Suppose p > 3/2. Then for each 7? G Wp~ 3/p (R n ) and m G Wp~ 3/p (W l ) 
there exists an extension n G Ey(l, 1) satisfying i](0) = 770, dtn(0) = n\ and 
the estimate 

IMIk»(i,i) < C (ll%|| w 4-3/p (Rn) + llJhll^-s/p^)) ■ 

(mj Suppose p > 3. T/ien for each rj € W / p 4 " 3/p (R") and 771 £ iy p 2 ~ 6/p (R n ) 
£/iere exists an extension 77 g E^(0, 1) satisfying 77(0) = 770, 9*77(0) = 771 and 
i/ie estimate 

IMIe*(o,i) < c (\\Vo\\ w ^-s/p {Mn) + IMI^-6/*> <■»„•,) ■ 
Proof, (i) Let 1 < p < 00. We claim that 

nit) := (2c-*( 1 - A ^ - c- 2 ^ 1 -^)^ (77) 
satisfies the properties asserted in (i). We have 

for fc = 1,2. It is a consequence of the mixed derivative theorem that the latter 
space is continuously embedded in Wp~ 1 ^ 2p ( J, Wp (R™)). This implies that 
^-fctci-A.)^ = „ Ax - )e -kt(i-A x ) m G W p 1_1 / 2p (J, Zp(R n )). 
Consequently, 

??e W ^-l/2p(J )Lp (R")) 

and we have that 

||77|| m/ 2-1/2 P(JLij(r „ )) < C|l 7 ?ll W - ? l(J iW -2-l/p (K „ ))nip(JJV 4-l/p (R „ :)) - 

The maximal regularity of (1 — A x ) on Wp~ 1 ^ p (W n ) and the embedding 

E|(l, 1) VK p 1 (J,VF p 2 " 1/p (M")) 

then yields 

IMIe*(1,1) < C||77o|| w 4-3/p (R „ ) . 

Obviously 77(0) = 7/0. If p > 3, the time trace of dtn is well defined and we also have 

9*77(0) = 0. This proves (i). 

(ii) Now suppose p > 3/2. Here we first set 

fj(t) := (e^ 1 ^) - e- 2t ( 1 - A ^)(l - A x )-%, (78) 
Then for 771 G Wp~ 3/p (R n ) we have that 

e-^-^il - A x )-% G % 1 (J,I¥ p 2 - 1/p (IR"))nL p (J,W p 4 - 1/p (K")) 
for fc = 1,2. By virtue of the embedding 

WpV, % 2 - 1/p (lR n )) n Lp{J, Wp 4 - 1/p (M")) M^- 1/2p (J, Wp 2 (R")) 
we obtain 

9ie - fet (l-A K ) (1 _ A^-l^ g Wl-l/2 P(J]Lp(M n ))) 
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hence that 

e -kt(i-A*) {1 _ Ax) -i T)i G W 2 -^ 2p (J,L p (R n )). 
By the same arguments as in (i) we obtain the estimate 

11*71^(1,1) < C !l r ?i||^-3/P (R „ r 
If fj denotes the extension constructed in (i) , then 

i] := fj + fj 

satisfies the regularity assertions in (ii). That 7/(0) = i]q and dtr)(0) = r\\ is obvious, 
(iii) Now we set 

fj(t) := (e-'( 1 - A ^ 2 - e-^-^Xl - A x )~ 2 r]i. (79) 
We have to check that c -kt(i-A x ) 2 ^ _ A^)- 2 ^ G E 2 (0,1). In view of r)i G 
Wp~ 6/p (R n ) we have that 

e -fct(i-A.) 2 (1 _ Ax )-2 r?i e W 2-2/ P ( R n^ n W e-VP^y 

From the embedding 

Wp{J, w 2 - 2/p (R n )) n L P {J, w£- 2/p (R n )) c — > W^ /2 - 1/2p (J, W^(R n )) 

we infer 

9 te -fe*(i-A,f (1 _ a,)- 2 ^ g V^ 2 - 1 / 2 ^^™)), 
and therefore that 

g-wfi-A.)^ _ Ax yz m € W* /2 - 1/2p (J,L p (R n )). 

Then rj := fj + fj satisfies all the assertions claimed in (iii), where fj denotes again 
the extension obtained in (i). □ 

Lemma 3.1 in conbination with [10, Lemma 6.4] yields the following result which 
provides a simultaneous extension for different regularity assumptions on the traces. 

Lemma 3.2. Let 3 < p < oo. T G (0,oo], and J = (0,T). For r/o and r/i there 
exists an (simultaneous) extension function rj such that 77(0) = r/o, dti](0) = r)i, and 



IMInsf.Co.o) < C (ll 7 ?oll w ,-2-2/p (R „ ) + IMI w i-3/ P(K 
if(Vo,Vi) G W 2 - 2/p (R n ) x W^ 3/p (R n ), 

IMIe»(o,i) < c (ll»«)llv^- s /»(R») + IMIw*- 6 /»(r»)) ' 
if(m,m) G Wp- 3/p (^ n ) x Wp' 6/p (R n ), and 

if {Vo,Vl) G Wp~ 3/p {R n ) x Wp~ 3/p (R n ), with C > independent 0/770 and 771. 

Proof. The idea for obtaining a simultaneous extension function as stated in the 
lemma is to employ a combination of the extension operators we used in Lemma 3.1. 
More precisely, we claim that 



n{t) := ^e^ 1 -^) 172 - e -^{i-A x )^ (V^ 1 "^) - e^ 1 "^) 
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satisfies all the properties asserted. Observe that (e _ ^*( 1_A<e ' )t>o is a bounded 
Co-semigroup and (1 — A^) - " is a bounded operator on WI(M. n ) for all r, a,/3 > 
and 1 < p < oo. Hence, 



e 



-W-^^l-AJ-" G J?(W r '(J,Wr(R n ))) 



for all s,r,a,/3 > and 1 < p < oo. If (770,771) G Wp 2/p (R") x Wp 3/p (R"), we 
therfore may estimate 

IMk(0,0) < c(||(2e-*( 1 -^) 1/2 -e- 2 *( 1 -^) 1/2 ) % ||^ (0i0) 

+ l|e- t(1 - A ^r/i|| E |(o,o)). 

By [10, Lemma 6.4] the remaining extension operators are known to lift the traces 
into the class E|,(0, 0), which implies 

IMIe*(o,o) < C (ll%|| w =-2/f (R „) + \\Vi\\ w ^/p {Rn) ) ■ 

Hence the first estimate is proved. If (r]o,T]i) G Wp~ 3 ^ p (W n ) x Wp~ 6 (M. n ), we 
interchange the roles of the semigroups in the definition of 77. In fact here we obtain 
as in Lemma 3.1 (hi), 

IMkcci) < cO^e-^-^-e-^-^MIlKo.l) 

+ ||(e-*(i- A *) 2 - e- 2 *< 1 ^) 2 )(l - A x )- 2 m \\E 2 T (0, 1)) 

< C + \\Vl\\ w 2-e/P {Kn) ) ■ 

Analogously we proceed in the third case. Here we treat the terms of type e _ ^*( 1_A ^) 1/2 
in front of 770 and the terms e -,3t ( 1-Ax ) and (1 — A^) -1 in front of 771 as bounded op- 
erators and gain the desired regularity by the remaining operators as in Lemma 3.1 (ii). 
A straight forward calculation also shows that 77(0) = 7/0 an d ^77(0) = □ 
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